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CeUs often perform computations in order to respond to environmental cues. A simple example 
is the classic problem, first considered by Berg and Purcell, of determining the concentration of 
a chemical ligand in the surrounding media. On general theoretical grounds, it is expected that 
such computations require cells to consume energy. In particular, Landauer's principle states that 
energy must be consumed in order to erase the memory of past observations. Here, we explicitly 
calculate the energetic cost of steady-state computation of ligand concentration for a simple two- 
component cellular network that implements a noisy version of the Berg-Purcell strategy. We 
show that learning about external concentrations necessitates the breaking of detailed balance and 
consumption of energy, with greater learning requiring more energy. Our calculations suggest that 
the energetic costs of cellular computation may be an important constraint on networks designed 
to function in resource poor environments, such as the spore germination networks of bacteria. 
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The relationship between information and thermody- 
namics remains an active area of research despite decades 
of study [l|-|4|. An important implication of the recent 
experimental confirmation of Landauer's principle, which 
relates the erasure of information to thermodynamic irre- 
versibility, is that any irreversible computing device must 
consume energy [2|, yj. The generality of Landauer's ar- 
gument suggests that this is true regardless of how the 
computation is implemented. A particularly interesting 
class of examples relevant to systems biology and bio- 
physics are intracellular biochemical networks that com- 
pute information about the external environment. These 
biochemical networks are ubiquitous in biology, ranging 
from the quorum-sensing and chemotaxis networks in 
single-cell organisms to networks that detect hormones 
and other signaling factors in higher organisms. 

A fundamental question is the relationship between 
the information processing capabilities of these biochem- 
ical networks and their energetic costs [5|, |6| . Energetic 
costs place important constraints on the design of phys- 
ical computing devices as well as on neural computing 
architectures in the brain and retina [7H9|, suggesting 
that these constraints may also influence the design of 
cellular computing networks. 

The best studied example of a cellular computation is 
the estimation of a steady-state concentration of chemi- 
cal ligand in the surrounding environment |10l4l2| . This 
problem was first considered in the seminal paper by Berg 
and Purcell who showed that the information a cell learns 
about its environment is limited by stochastic fluctua- 
tions in the occupancy of the receptors that detect the 
ligand [10|. In particular, they considered the case of a 
cellular receptor that binds ligands at a concentration- 
dependent rate k'^^ and unbinds particles at a rate fc|* 
(see Fig. [T]) and argued that cells could estimate chemical 
concentrations by calculating the average time a recep- 
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FIG. 1. A cellular network for the computation of 
an external ligand concentration. External ligands are 
detected by a receptor that can exist in two conformations: a 
high-activity "on" state and a low-activity "off". Receptors 
switch between states at rate fcoff and fcon- Receptors in state 
s=on,off can post-translationally modify (i.e. phosphorylate) 
a downstream protein at a rate fcl. This modification is lost 
(i.e. dephosphorylated) at a rate fei . 



tor is bound during a measurement time T ^ 1. Re- 
cently, however, it was shown that the optimal strategy 
for a cell is instead to calculate the average duration of 
the unbound intervals during T or, equivalently, the to- 
tal time that the receptor was unbound during T. This 
later computation implements Maximum Likelihood Es- 
timation (MLE) [12|. In these studies, the biochemical 
networks downstream of the receptors that perform the 
desired computations were largely ignored because the 
authors were interested in calculating fundamental lim- 
its on how well cells can compute external concentra- 
tions. However, calculating energetic costs requires us to 
explicitly model the downstream biochemical networks 
that implement these computations [13| . 



Here, we consider a simple two-component biochemical 
network that encodes information about ligand concen- 
tration in the steady-state concentration of the activated 
form of a downstream protein (shown in Figure [T]). Such 
two-component networks are common in bacteria and of- 
ten used to sense external signals with receptors phos- 
phorylating a downstream response regulator IJ]. Re- 
ceptors exist in an active 'on' state and an inactive 'off' 
state. For simplicity, as in previous works |10l - ll2l | , we as- 
sume that the binding affinity in the on state is extremely 
high so that the ligand-bound receptors are always in the 
on state and unbound receptors in the off state. Recep- 
tors can switch between the off state and on state at a 
concentration- dependent rate k^^ and from the on state 
to the off state at a concentration-independent rate fc™ . 
Receptors convert a downstream protein from an inactive 
form X to an active form X* at a state-dependent rate 
fc|, with s = on, off. The proteins are inactivated at a 
state- independent rate fci. 

The phosphorylation rates in the off state are small 
but must be nonzero for thermodynamic consistency [15| . 
^2 includes non-specific phosphorylation due to other ki- 
nases as well as contributions from the reverse reactions 
of the phosphotases. The inactivation rate sets the ef- 
fective measurement time T oa k^ since it is the rate 
at which information encoded in downstream proteins is 
lost due to inactivation. In order to compute external 
concentrations accurately, the measurement time must 
be much longer than the typical switching times between 
receptor states, fci ^ fc™,^;^^. We show below that this 
simple network actually implements a noisy version of 
the original Berg-Purcell calculation and discuss the re- 
lationship between information and power consumption 
in the context of this network. 



STEADY-STATE PROPERTIES 



The variance of the Langevin terms is given by the Pois- 
son noise in each of the reactions 

{lln{t)Vn{t')) = (fc^Pon + k"^^ ^ {I - p,^) + k,n)S{t - t') 
{Vr(t)Vr(t')) = (fcf(l - Pon) + fcrPon)<5(i - t'), (3) 

with 6{t — t') denoting the Dirac-delta function and 
barred quantities denoting the mean steady-state values 

[3 [13 

At steady-state, we can calculate the mean probability 
and number of proteins by setting the time derivative in 
Eq. [2] equal to zero and ignoring noise terms, yielding 



1 - Poff = 



/.-off 



K2" + K2' 



and 



(i^r - ^2")Po„ + k: 



-off 



(4) 
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where we have defined the dimensionless parameters 
K!- — kj/ki with j — 2,k and s = on, off. For the bi- 
ologically realistic case K2^ <^ i^^Pon, as expected, the 
mean number of proteins is proportional to the kinase 
activity in the on state times the probability of being in 
the on state, n « i^2"Pon- One can also calculate the 
variance in protein numbers (see Appendix) 



(iSn)^)=n+{AKrf 
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The first term on the right hand side of the equation 
arises from Poisson noise in the synthesis and degradation 
of activated protein, whereas the second term is due to 
stochastic fluctuations in the state of the receptors. 

In addition to the mean and variance, we will need 
the full, steady-state probability distributions for n to 
calculate power consumption. The steady-state distribu- 
tion can be calculated from the master equation for the 
probability, Ps{n), of having n active proteins with the 
receptor in a state s. 



The deterministic dynamics of the biochemical net- 
work in Figure [T] is captured by simple rate-equations 
for the mean number of activated proteins. We can aug- 
ment these equations to account for stochastic fluctua- 
tions within the Linear-noise approximation by adding 
appropriate Langevin noise terms. In what follows, we 
assume that proteins are abundant and ignore saturation 
effects. The dynamics of the cellular circuit is described 
by a pair of Langevin equations for the probabilities pon 
and Poff of the receptor to be in the on and off state, 
respectively, and the number of activated proteins n. 
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kf{l~Pon)-krPon + vAt) (1) 

fc^Pon + k°/^\l - pon) - km + ?7„(t). (2) 



dt 



ki{n + l)ps{n + 1) + k^psin - 1) 

+ klps{n)~{kin + k^ + kl)ps{n) (7) 



with s = off (on) when s = on (off). At steady-state, the 
left-hand side of Eq. [7] is zero and 

Klpsin) = -{n + l)ps{n + 1) - K^Psin - 1) 

+ {n + K^ + Kl)ps{n). (8) 

This equation is similar to those found in [18|, |l9[ and 
can be solved via a generating function approach. Define 
a pair of generating functions. 



Gs{n) = J2Ps{n)z\ 



(9) 



n=0 



with s — on, off. We can rewrite 
generating functions as 



in terms of the 



[{z - l)d, - K^{z - 1) + Kl] Gsiz) = KlGsiz). (10) 

This equation must be supplemented by initial conditions 
for the Gs{z). These follow from the observation that 
Gon(l) = Pon and G'off(l) = PoS = 1 -Pon with pon given 
by Eq. U) As shown in the Appendix, this equations can 
be solved exactly and yields 
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Kie 
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(11) 
where AKI = AK2 — AK2 and iFi{a; 6; z) is the conflu- 
ent hypergeometric function of the first kind. As a check 
on this expression, we can compute the variance of n di- 
rectly from Eq. [TT] and it is in agreement with ^ (see 
Appendix). 

Depending on the parameters, the steady-state dis- 
tributions have two qualitatively distinct behaviors (see 
Fig. [2]). In the slow switching regime with fc§^ <^ fc™ 
and A;™, kf^ <^ fci, receptors switch at rates much slower 
than the protein deactivation rate fci. This gives rise to 
a bimodal distribution of activated proteins that can be 
roughly thought of as a superposition of the probability 
distribution of activated proteins when the receptor is the 
on and off state. As ^2" approaches ^2^, the distributions 
in the two states merge and the overall probability distri- 
bution becomes unimodal. On the other hand, in the fast 
switching regime, where fc°",fc°* <^ fci, the distribution 
of activated proteins is always unimodal. In this regime, 
the measurement time T oc fcj" is much longer than the 
average time a receptor is in the on or off state, and 
the biochemical network 'time-averages' out the stochas- 
tic fluctuations in receptor states. In what follows, we 
restrict our considerations to this latter regime. 



QUANTIFYING LEARNING 

The biochemical circuit in Figure[T] "computes" the ex- 
ternal concentration of a chemical ligand by implement- 
ing a noisy version of Maximum Likelihood Estimation. 
As emphasized by Berg and Purcell in their seminal pa- 
per [10| , the chief obstacle in determining concentration 
is the stochastic fluctuations in the state of the ligand 
binding receptors. Berg and Purcell argued that a good 
measure of how much cells learn is the uncertainty in ex- 
ternal concentration as measured by the variance of the 
estimated concentration (Sc)'^ . Berg and Purcell assumed 
that the cell computed the average receptor occupancy by 
time-averaging over a measurement time T, and showed 
that E^ 
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FIG. 2. Top. Slow switching regime, fci > k'r,kf. To- 
tal probability (black solid line), probability when receptor 
is in the on state (blue dash-dot line), probability when re- 
ceptor is in the off state (red dashed line). Middle. Fast 
switching regime, fci <^ k'^,k1 . Total probability (black 
solid line), probability when receptor is in the on state (blue 
dash-dot line), probability when receptor is in the off state 
(red dashed line). Bottom. The uncertainty in ligand con- 
centration, (Scrms/c) as a function of fci with mean number 
of active proteins n = 25 (dashed red line) and n = 100. This 
can be compared to the Berg-Purcell result (solid black line). 
Parameters: kf = 0.01, fc?", fc^ = kf = 1. 



with k1^ — k^c, fc™ = fc- and independent of c, and 



Nh the number of binding events during the time T. It 
was later shown that cells could compute concentration 
more accurately by implementing Maximum Likelihood 
Estimation (MLE) with [H, ^ 



{ScMhf 



1 JSCBP)' 

2"" C2 



(13) 



The decreased uncertainty is a results from the fact that 
MLE ignores noise due to unbinding of ligands from the 
cell. 

To quantify learning in our biochemical circuit, we fol- 
low Berg and Purcell and estimate the fluctuations in 
((5c) 2 as 



(Sc)' 



e|) iSnr 



(14) 
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with (Sn)'^ = (n^) — n^. Substituting k^ = k+c and 
fcf" = fc- and computing the derivative using Eq. [S] 



gives 



dn 
dc 



{PonPoff^K2) 



Substituting Eq. [5] and IH into Eq. [TH yields 



L,0ff 



> ki (equivalently K^^ > K^*^ > 1), and the av- 
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erage number of activated proteins is large, h ^ ^2 ^ 
Qr\ 1, the expression above reduces to ((5c)^/c^ w 2/iV{, in 
agreement with the Eq. [121 
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(16) 
Similar to the linear noise calculation, the first term on 
the right-hand side arises from the Poisson fluctuations in 
activated protein number whereas the second term comes 
from the stochastic fluctuations in the state of receptors. 
Figure [2] shows the uncertainty, (5c) ^/c^, as a function of 
the degradation rate of activated protein, fci when n = 25 
and n = 100 and A:§" > kf . 

By identifying the degradation rate with the inverse 
measurement time, fci = 2r~^, we can also compare the 
results with Berg Purcell. The factor of two is due to the 
slight difference in how the variance of the average re- 
ceptor occupancy is calculated for a biochemical network 
when compared to the original Berg Purcell calculation 
[20| . As shown in Fig. [51 when fi is increased, the Poisson 
noise in protein production is suppressed and the perfor- 
mance of the cellular network approaches that due to 
Berg-Purcell. To make the connection with Berg-Purcell 
more explicit, it is helpful to rewrite Eq. [161 in terms of 
the average number of binding events, N},, during T 



POWER CONSUMPTION AND ENTROPY 
PRODUCTION 



We now compute the energy consumed by the circuit 
in Figure [1] as a function of the kinetic parameters. To 
do so, we exploit the fact that dynamics of the circuit 
can be thought of as a nonequilbrium Markov process 
(see Fig. [3]). A nonequilibrium steady-state (NESS) nec- 
essarily implies the breaking of detailed balance in the 
underlying Markovian dynamics and therefore a non-zero 
entropy production rate. The entropy production rate is 
precisely the amount of power consumed by the biochem- 
ical circuit to maintain the nonequilibrium steady-state. 
Thus, by calculating the entropy production rate as func- 
tion of kinetic parameters, we can calculate the power 
consumed by the biochemical network implementing the 
computation. 

Consider a general Markov process with states labeled 
by a and transition probability from a to a' given by 
k{(T,a'). Defining the steady-state probability of being 
in state a by Pa-, the entropy production rate for a NESS 
is given by [21], 
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When the measurement time is much longer than 
the timescale for fluctuations in the receptor number, 

I 



^'^ S:^ uf \i t 'M k{a,a') 
-=^P(a)fc(a,a)log^^-;-^, 



(18) 



For the biochemical network described by Eq. [3 this 
becomes. 
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Since the receptors are in thermodynamic equilibrium, 
from detailed balance we know that 



/C/i 



^p.(n)fc|log-| = 0. 



(20) 



so that 




s=on,off n ^ 


- n log 

n 



(21) 
where K^ — fc|/fci. The steady-state distributions Ps{n) 
follow from Eq. [TTl The physical content of this ex- 
pression is summarized in Fig. [31 The expression states 



that any non-zero cyclic flux must necessarily produce 
entropy. If it didn't, one would have a chemical version 
of a perpetual motion machine. Figures [31 and 01 show the 
power consumption as a function of AK2 = ^^2" ^ ^2^ 
and fci . Notice that the power consumption tends to zero 
as both these parameters go to zero. We cannot, how- 
ever, set fci = identically because there then no longer 
exists a steady-state distribution. 




FIG. 3. Top. The probabilistic Markov process underlying 
the circuit in Fig. [l] Any non-zero cyclic flux (depicted in 
red) results in entropy production and power consumption. 
Bottom. Power consumption (solid black line) and uncer- 
tainty (dashed purple line) as a function of Afc2 — fc™ — fc2 
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FIG. 4. Total energy per independent measurement (Power 
xfcj"^) as a function of fei when n = 25. (Inset) Power as a 
function of fci . 
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probabilities of being in the on or off state. 

Second, as shown in Fig. |4l the power consumption 
of the circuit tends to zero as fci — > 0. This is consis- 
tent with Landauer's principle: entropy production stems 
from erasing memory in a computing device. The num- 
ber of activated proteins serves as a memory of ligand 
concentration which is erased at a rate k\. Thus, as the 
erasure rate of the memory tends to zero, the device con- 
sumes less energy, as expected. Yet despite the fact that 
the power consumption tends to zero as fci decreases, 
the total energy consumed per measurement, namely the 
power times the measurement time, T ~ 2fcj~ , still in- 
creases (see Fig. 14]). Thus, learning more requires con- 
suming more total energy despite the fact that power 
consumption is decreasing. In effect, one is approaching 
the reversible computing limit where memory is erased 
very infrequently. Note, however, that when erasure is 
performed infinitely slowly, fci = 0, the system no longer 
has a NESS and our formalism does not apply. 

Finally, we note that one of the important open prob- 
lems in our understanding are the constraints placed on 
the measurement time T . In principle, cells can always 
learn more by measuring the environment for longer pe- 
riods of time. However, in practice, these measurement 
times tend to be quite short. There are a number of 
constraints that can set this measurement time including 
rotational diffusion [lO[ and the restrictions placed on 
motility. Here, we highlight another restriction that may 
be important in resource-starved environments: sensing 
external concentration necessarily requires cells to con- 
sume energy. 



DISCUSSION AND CONCLUSION 



ENERGETICS, INFORMATION, AND 
LANDAUER'S PRINCIPLE 



We now highlight the fundamental connection between 
the energy consumed by the network and the information 
the network learns about the environment, and briefly 
discuss its relation to Landauer's principle. First, note 
that learning information about the environment requires 
energy consumption by the network. This can be seen in 
Fig. [3] which shows that as Afc2 -^ 0, the uncertainty 
about the concentration tends to infinity. This can be 
made concrete by noting that the entropy production 
(Eq. W\\ is zero if and only if l\ki — (see Appendix) . In 
conjunction with Eq. 1161 which diverges as Afc2 ~^ 0, this 
implies that learning requires consuming energy. Phys- 
ically, in the limit where Afc2 — 0, the dynamics of the 
Markov process in Fig. [3] become "one-dimensional" and 
the dynamics obey detailed balance. At the same time, in 
this limit, the number of downstream proteins becomes 
insensitive to external ligand concentrations, since all in- 
formation about concentration is contained in the relative 



Cells often perform computations using elaborate bio- 
chemical networks that respond to environmental cues. 
One of the most common simple networks found in bac- 
teria are two-component networks where a receptor phos- 



phorylates a downstream response regulator IJ] . In this 
work, we have shown that these simple two-component 
networks can implement a noisy version of the Berg- 
Purcell strategy to compute the concentration of external 
ligands. Furthermore, by mapping the dynamics of the 
biochemical network to Nonequilibrium Steady-States in 
Markov processes, we explicitly derived expressions for 
the power consumed by the network and showed that 
learning requires energy consumption. Taken together, 
these calculations suggest that, much like man-made and 
neural computing ^, ,7,-1^] , energetic considerations may 
place important constraints on the design of biochemical 
networks that implement cellular computations. They 
also suggest a fundamental relationship between the effi- 
ciency of cellular computing and the energy consumption. 
Bacterial cells such as Bacillus suhtilis can sporu- 
late during times of environmental stress and remain 



metabolically dormant for many years. While sporula- 
tion is relatively well understood, the reverse process of 
germination is much more difhcult to study. One cur- 
rent model for how a spore knows when to germinate 
in response to external cues involves integrating the sig- 
nal and triggering commitment when an accumulation 
threshold is reached 22|, |23|. This corresponds to the 
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limit of vanishingly small fci in our model, so that power 
consumption is minimized at the expense of retaining the 
entire integrated signal. Our results indicate that this be- 
havior may be due to the extreme energetic constraints 
imposed on a metabolically dormant spore, rather than 
an evolutionarily optimized strategy. 

An important insight of this work is that even a sim- 
ple Berg-Purcell strategy for sensing external concentra- 
tions requires the consumption of energy. It is likely that 
more complicated strategies that increase how much cells 
learn, such as Maximum Likelihood, require additional 
energetic inputs. For example, it was argued in [20| that 
MLE can be implemented by a network similar to the 
perfect adaptation network where bursts are produced in 
response to binding events. These bursts break detailed 
balance and therefore require energy consumption. It 
will be interesting to investigate further how the trade- 
off between learning and energy consumption manifests 
itself in the design of computational strategies employed 
by cells. 

In this work, we restricted ourselves to the simple case 
where cells calculate the steady-state concentration of an 
external signal. In the future, it will be useful to general- 
ize this to other computations such as resp onding to tem- 
poral ramps [20| and spatial gradients 2J,|25|. It will also 
be interesting to understand how to generalize the con- 
siderations here to arbitrary biochemical networks. An 
important restriction on our work is that we reduced our 
considerations to nonequilibrium steady-states. It will 
be interesting to ask how to generalize the work here to 
biochemical networks with a strong temporal component. 
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Appendix: Variance from the Linear-noise Approximation 

Consider the Langevin equations for the ODE's 

fcf(l-Po„)-fc4>o„ + ^r(i) (22) 



dPon _ , oS/ 



dt 

5 = kTPon + kf{l ~ Pon) - hn + r;„(i) (23) 

at 



We can linearize (with bar denoting average) to get 

-(fcf + kDSpon - Tp'Spon + V^W (24) 



d^Pon ,,o« 



where 



dt 

-^ = (fc°" - kf)Spon -ki5n + rjn (25) 



{'np{t)r,p{t)) - fcf (1 - Pon) + k^PoU = 2ikf + fcr)Po„(l - Pon) = 2Tplpon(l " Pon) (26) 



and ('r]n(t)'>]n{t)) = 2kin. Now we Fourier transform the above equations and use the fact that 

{{5nf)= f'^{6h{uj)5h*{^)) (27) 



to get 



((^")^) - ^ + ^^f^'. 'l-i - ^ + i^Kr) \J;±'.„ (28) 



where we have used K^ — kf/ki. 

Appendix: Generating function for probability distribution 

We start with the equation for the generating functions in the main text, 

[{z - 1)5, - if|(z - 1) + Kl] Gs{z) = KlGs{z)- (29) 

Adding the equations for s =on, off and dividing through by [z — 1) gives 

(a, - Ar)Gon(^) = -{d,- Kf)Goff{z). (30) 

To proceed further, it is useful to define the quantity Hs{z) related to the generating functions Gs{z) by 

G,{z) = e''i'Hs{z). (31) 



It is clear that 



Thus, we can rewrite (|30| as 



(9, - K'^)Gs{z) = e^'i'd^H^iz) (32) 



e^S"^a,i/on(^) = -e^^° 'd^Hosiz). (33) 

Similarly, we can rewrite p9p as 

(z - l)e''i'd,H,{z) + Kle^'i'H.iz) = Kle''^'-H,{z) (34) 



Multiplying the equation by e ^2^, taking the derivative withe respect to z, substituting ([33]), and defining AK2 — 
one has 

d,H,{z) - AK^iz - l)d,Hs{z) + (z - l)dlH,{z) - K'^AK^Hs{z) + Kid,H,{z) = ^Kld,H,{z) (35) 



K2 — K2S one has 
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Regrouping terms one has 

(z - l)d^,Hs{z) + (1 - AK^{z - 1) + Kl + Kl)d,H,{z) - KlAK^H,{z) = 0. (36) 

Defining 

u = AK^{z~l), (37) 

one can rewrite (I36p in terms of u as 

udlHs{u) + (1 + if| + Kl - u)duHs{z) - KlHs{u) = 0. (38) 

This is just the confluent hypergeometric equation. We can immediately write the solutions in terms of confluent 
hypergeometric functions of the first and second kind. In particular, the general solution to this equation is given by 
(assuming that we need a power series in u with integer cofficients) the confluent geometric function: 

Hs {u) ^c.iFi {Kl- 1 + KI + Kl; u), (39) 

with Cs a constant of integration. Thus, we have 

G,(z) = c,e^^^^Fi(if|; 1 + if | + X|; AX|(z - 1)). (40) 

To determine the constants, notice that 

GAl) = CsC^'^ ^Ps = — ^-- (41) 



so that 



This gives the final expression in the main text 



e-^^Kl 



Gsiz) = ^^l ^ ^,' iF.iKl; 1 + KI + Kl; Aif|(z - 1)). (43) 

Appendix: Variance from full probability distribution 

We can calculate the variance (Sn)'^ directly from the generating functional. This serves as a check on our formalism 
and our expressions. In terms of the generating functions, 

(Sn)^ - Y,idzzd,G,{z))U^, -n'= J2i^d^G,{z))U^, + n - n' (44) 

s s 

where we have used 

zd^ iFi{a,h; z) = z- iFi(a+, 6+; z). (45) 

a 



Taking the second derivative and setting z = 1 yields 



{zdlGM)l=-- ^' 



Kl + Kl 
After some algebra, one has 



, AK^.Kl ( , Ai^l 



l + Kl + Ki\ ' 2 + Kl + Ki 



(46) 



LTon I LToS 

Y,{zdlGs{z))U=i ^ Pon{Kn'+Poff{Kf)' - ponPoffiAKir , ^ 'Z^ toff (47) 

This gives 

[8nf = n + p.^{KTf + p.^{Kff - p„„p„ff(AX|)2--i^-l-i^ - n^ (48) 

= n + p.AKTf + Po^{Kff - ponPofF(AJfl)^ ^ Z^ toff - {VonKT-r Po^Kff (49) 



Appendix: Energy consumption is required for signaling 

We show is that ^ = if and only if AK2 = 0. To do so we wiU use the expressions for the probabihty distribution 
PTjl and average entropy production (I^T]) . Let us start by writing (pi]) as 

/cfi— = - ((n)o„ logi^r + (n)off logiff ) + Y, (K^Ponin) logi^r + Kfposin) \ogKf) (51) 

n 

+ J2 ([Ponin + 1) - K°y,^{n)] log (n + 1)) + J] {[Pos{n + 1) - i^fpoffCn)] log (n + 1)) , (52) 

n n 

where {n)s = ^^^ nps{n) for s = on, off. We will now examine each of these terms one by one. Notice that 

(n}.=9.G.(z = l) = -^^2^, (53) 

so we have that 

{{nUlogK^ + (nUhgKf) =J2^^^\ogK^. (54) 

S ^ ■^ 

Notice that since X!-Ps(^) — *^s(l)i ^^ second term in the sum above is just 

^ {KTPoAn) logi^r + ^2°Vff(^) ^ogKf) = Y. ks^Tk-s ^°S^2- (55) 



Kl + XI 



Substituting these expressions into ([57)) yields 



/ = - E ]^r^ log^2 + E ([P°n(" + 1) - ^2>o„(n)] log (n + 1)) (56) 

s ^ "^ n 

+ Yl ([P°ff(" + 1) - ^^'2 Vff(")] log {-n + 1)) • (57) 

n 

In order for this to equal zero, we know that 

p,(n+l) = -^p,(n), (58) 

n + 1 

since each of the terms in the sum above must equal zero individually. This implies that the generating functional in 
detailed balance, G'f^(z), must be given by 

G^^ ex e^2°^ (59) 

Comparing with pT|) implies that IS.K2 = 0. In addition, the first term of ([57]) also disappears in this case. Thus, we 
have detailed balance if and only if A_ftr| — 0. 



